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ABSTRACT 

We calculate the running of the three coupling constants in cosmological, topo- 
logically massive 3d gravity. We find that u, the dimensionless coefficient of the 
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stant and Newton's constant depends on v, and for any positive v there exist both a 
trivial and a nontrivial fixed point. 
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1 Introduction 



Einstein gravity in three dimensions with or without cosmological constant and sup- 
plemented with a gravitational Chern-Simons (CS) term has been a subject of many 
investigations since long pQ. This theory is of considerable interest, as the presence 
of the cosmological term makes it possible to have a black hole solution, while the 
Chern-Simons term is responsible for the presence of a single propagating massive 
graviton. For the latter reason, this model is referred to as the topologically massive 
gravity (TMG). 

The properties of this theory have been the subject of intense scrutiny. For a 
generic value of the CS coupling, either black hole states (if G < 0) or graviton states 
(if G > 0) will have negative mass. In [2] it was observed that in the case G > 0, if 
one chooses a critical value of the CS coupling, called the chiral point, the negative 
mass graviton mode can be confined to propagate only on the boundary, provided 
that suitable boundary conditions are imposed. The prospects of AdS/CFT duality 
are enhanced in this framework. There have also been alternative views in which the 
bulk graviton is maintained but the negative energy black hole solution is viewed as 
being possibly irrelevant [3]. 

TMG appears to be renormalizable [1H6]. However, not much is known about 
its renormalization group (RG). For example, does the chiral point have any special 
properties? Another motivation for addressing this issue comes from the asymptotic 
safety approach to quantum gravity [7JIE] , which requires the existence of a fixed point 
(FP) with finitely many UV attractive directions. Such a nontrivial FP has been found 
in the Einstein-Hilbert truncation of gravity in any dimension d > 2 [9HTT]. but in 
d = 3 this theory does not have propagating degrees of freedom. A nontrivial FP in 
TMG could provide an example of a three dimensional asymptotically safe theory of 
gravity with physical degrees of freedom. 

Our purpose in this paper is to calculate the one loop beta functions for TMG, to 
determine the RG flow of the couplings and to discuss their FPs. We find that many 
subtleties and technically nontrivial issues arise due to the presence of the Lorentz 
Chern-Simons. For example, the time honored heat kernel expansion method is not 
well suited in this case, essentially due to the fact that the eigenvalues of the cubic 
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in derivative wave operator present in the theory is a cubic polynomial in the lowest 
weight labels, as opposed to being quadratic. There is the delicate matter of how 
to choose the cut-off scheme to do these computations as well. Our results are best 
discussed after spelling out all these issues, and are summarized in the Section 7. As 
we shall see, one of the key results we will obtain is that the dimensionless coefficient, 
is, of the Chern-Simons term has vanishing beta function. Moreover, we find that the 
flow of the cosmological constant and Newton's constant depends on the values of v, 
and for any positive v there exist both a trivial and a nontrivial FP. Various aspects 
of these results are also discussed in Section 7. 

This paper is organized as follows. In Section 2, we review briefly the Wilsonian 
approach to the study of Renormalization Group (RG) equations. In Section 3, we 
determine the gauge fixed inverse propagator of TMG on a maximally symmetric 
background. In Section 4, we define the IR cutoff of each spin component separately 
and we give the harmonic expansion of the fields. In Section 5, we evaluate the 
functional traces that arise in the computation of the one-loop beta functions, using 
the Euler-Maclaurin formula. The resulting beta functions and their FPs are given 
in Section 6. In Section 7, we first summarize our results, and then discuss various 
aspects of them a number of subsections. Four appendices are provided, of which 
the first two contain useful lemmas and the details of the Euclideanization procedure 
and harmonic expansions on S 3 . In Appendix C, we address detail the role of the 
tadpoles in an off-shell beta function computation, and in Appendix D we give a heat 
kernel evaluation of the beta functions of pure gravity with cosmological constant, 
which provide a useful check of the techniques used in Section 5. 

2 Wilsonian Method For Computing the RGE 

There is no single approach to the study of Renormalization Group (RG) equations. 
Here, we shall adopt a Wilsonian method. The central lesson of Wilson's analysis of 
QFT is that the "effective" (as in "effective field theory" ) action describing physical 
phenomena at a momentum scale k can be thought of as the result of having integrated 
out all fluctuations of the field with momenta larger than k. Since k can be regarded 
as the lower limit of some functional integration, we will usually refer to it as the 
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infrared cutoff. The dependence of the "effective" action on k is the Wilsonian RG 
flowB 

There are several ways of implementing this idea in practice, resulting in several 
forms of the RG equation [13]. In the specific implementation that we shall use 
[T3] , instead of introducing a sharp cutoff in the functional integral, we suppress the 
contribution of the field modes with momenta lower than k. This is obtained by 
modifying the low momentum end of the propagator, and leaving all the interactions 
unaffected. We describe first the general idea, and comment on its application to 
gravity later. We start from a bare action S[<p] for some fields 0, and we add to 
it a suppression term ASk[4>] that is quadratic in the field. In flat space this term 
can be written simply in momentum space. In order to have a procedure that works 
in an arbitrary curved spacetime we choose a suitable differential operator O whose 
eigenfunctions <p n , defined by Oip n = \ n ip n , can be taken as a basis in the functional 
space we integrate over: 

00) = ^0n¥?nO) , (2.1) 
n 

where <p n are generalized "Fourier components" of the field. (We will use a notation 
that is suitable for an operator with a discrete spectrum.) Then, the additional term 
can be written in either of the following forms: 

= \ I dx<P{x)R k {0)<P{x) = ^Y,4>nMK). (2.2) 

n 

The kernel Rk{0) will also be called "the cutoff". It is arbitrary, except for the 
general requirements that Rk(z) should be a monotonically decreasing function both 
in z and k, that Rk{z) — > for z> k and Rk{z) ^ for z <C k. These conditions are 
enough to guarantee that the contribution to the functional integral of field modes 
(f) n corresponding to eigenvalues A n <C k 2 are suppressed, while the contribution of 
field modes corresponding to eigenvalues A n ^> k 2 are unaffected. We will further 
require that R k (z) — » k 2 for k — > 0. We define a fc-dependent generating functional 

2 At this general level of discussion, it is not necessary to specify the physical meaning of k: for 
each application of the theory one will have to identify the physically relevant variable acting as 
k. In scattering experiments k is usually identified with some external momentum. See |12l for a 
discussion of this choice in concrete applications to gravity. 
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of connected Green functions by 

-W k [J] 



e 



J £>0exp - AS k [(p] - J dxJ^ (2.3) 

and a modified /c-dependent Legendre transform 

T fc [0] = W k [J] - JdxJ<j>- AS k [cf)} , (2.4) 

where AS k [4>] has been subtracted. The functional T k is sometimes called the "effec- 
tive average action" , because it is closely related to the effective action for fields that 
have been averaged over volumes of order k~ d (d being the dimension of spacetime). 
The "classical fields" SW k /S J are denoted again for notational simplicity. In the 
limit k — > this functional tends to the usual effective action r [</>], the generating 
functional of one-particle irreducible Green functions. It is similar in spirit to the 
Wilsonian effective action, but differs from it in the details of the implementation. 

The average effective action used at tree level, should give an accurate 

description of processes occurring at momentum scales of order k. In the spirit of 
effective field theories, one assumes that T k exists and admits a derivative expansion 
of the form 

oo 

r.(^«) = EEft (n) W « w W , (2.5) 

n=0 i 

where gl n \k) are coupling constants and 0\ n ^ are all possible operators constructed 
with the field and n derivatives, which are compatible with the symmetries of the 
theory. The index % is used here to label different operators with the same number of 
derivatives. At one loop, T k is given by 

rl 1) = ^ + ^Trlog(^+i? fc ) (2.6) 

where denotes the second variation of the bare action. The only dependence on 
k is contained in the cutoff term R k and therefore 

kd -f + ^y 1 ^ ■ < 2 - 7 > 

Note that due to the properties of R k (z), the factor k^- goes to zero for z > k 2 . 
Thus, even though R k (z) is introduced in the functional integral as an infrared cutoff, 
in the evaluation of the beta functions it acts effectively as UV cutoff. One can extract 
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from this functional equation the one loop beta functions, by making an ansatz for 
Ik of the form (12.51) . possibly containing a finite number of terms, inserting it in the 
left hand side of (|2.7p . expanding the right hand side and comparing term by term. 

When applied to scalar, fermion and gauge field theories this method reproduces 
the well known one loop beta functions. In this paper we shall apply this technique to 
TMG. In general, in applying Wilsonian ideas to gravity one has to confront the fact 
that the definition of a cutoff generally makes use of a metric and since the metric is 
now to be treated as a dynamical field it is less clear what one means by cutoff. We 
will follow [15] and use the background field method, thereby effectively replacing the 
dynamical metric g^ u by a spin two field h^ u propagating in a fixed but unspecified 
background g^ u . In addition to appearing in the gauge fixing, the background metric 
can be used to unambiguously distinguish what is meant by long and short distances, 
and hence low and high energies. In practice, it allows us to write a cutoff term that 
is purely quadratic in the quantum field, as required in the preceding discussion. The 
effective average action is then a functional of two fields T^g^, h^ u ); in practice, aside 
from the gauge fixing and cutoff terms, which are quadratic in h^, we will only need 
the restriction of the functional Tk to the space where {h^J) = 0. This is sufficient to 
extract beta functions for the couplings we are interested in. 



3 The Quadratic Gauge Fixed Action for TMG 

Topological massive gravity is described the action 

s = z j <p x v=g [r-2K + 1 s^ri (d^i p + fr^ry ) , (3.1) 

where Z = []. At this point we need not specify the signs of the couplings G, A 
and fi. It will be useful to define the dimensionless combinations 

v = [iG \ r = AG 2 ; = n/^\K\ . (3.2) 

Note in particular that is the coefficient of the CS term. We will expand around 
background metric g^ v as g^ u = g^ v + h^ v . Next, we add to the action the gauge 
fixing term 

Sgf = -^J d 3 x^ X ^ u Xv , (3.3) 
3 We use the conventions R^ a = d^T^, - d v Y/ a + T/ T T„ T a - r^ r r M T ff and R„ v = R p / V . 
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where 



Xu = - ^d v h . (3.4) 



Note that — 1 corresponds to the familiar de Donder gauge. The ghost action 
corresponding to the gauge (13 .4p is given by 

S gh = - J Sx^-gC^ [b v p + ^y^V„ V + R^j C v . (3.5) 

The one loop average effective action 1^ contains terms of the same form as S + 
Sqf+ Sgh^except that all the couplings in S are now to be interpreted as renormalized 
couplings Q Their beta functions will be obtained from equation (I2.7p . which for this 
theory has the form 

k ^f _ I Tr ( S 2 (S + S GF ) \ - 1 dR^ ( 5 2 S gh „ V 1 dW\ 

dk 2 V Sh^Sh^ J dk \5C»5C V + V dk • 

(3.6) 

Expanding in powers of h^, and discarding total derivative terms, the quadratic 
part of the action is given by 

+ Sgf = \z j d 3 x^=g V ( D - y + 2A ) htiU + 2(1 a " a) VVV P r 

+ (2- hV»V u K v - ( 1 - ^±11!^) hUh + -h(R - 6A)h 

V a J ' V 8a J 6 

(3.7) 



The raising and lowering of indices is understood to be with the background metric, on 
which we drop the bar for notational simplicity. For our purposes it will be sufficient 
to consider maximally symmetric backgrounds, for which 

R R 

Rfii/pa ~^{.9^ip9u(j QiurQvp) 1 Rpu ~^9^f ■ (3-8) 

Here R = ±6/£ 2 , with the + sign for de Sitter, and the — sign for anti-de Sitter 
space, and £ is the "radius". It is important to note that, in computing the gauge 
fixed quadratic action and the beta functions, we do not use the equation of motion 
which for this class of metrics would imply the relation i = 1/ a/|A|. As we shall see 
later, this is necessary to compute the beta functions for A and G separately. 



We allow only the couplings in S to run, whereas the gauge fixing parameters will be kept fixed. 
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In order to achieve partial diagonalization of the inverse propagator we decompose 
the quantum fluctuation h^ u into irreducible parts: 

hfiu = fyL + + V,^ + V v V u a - -g^Ua + ^g^h , (3.9) 

where h = g^h^, £ M satisfies V a £a = and hZ v satisfies g^ u h^ u = and V A //£, = 0. 
It is important to note that in passing from a path integral over the field to one 
over the fields (h£ v , £, a, h), one has to take into account the appropriate Jacobian 
factors. These Jacobians will be exactly canceled by other Jacobians arising from 
suitable redefinitions of £ and a, as we will see later (see (14.2j) ). 

Inserting (13.91) in (13. 7p . and using the lemmas in Appendix A, we find 



S (2) + S GF = \z k J d"x^-g\hl v {u- 2 -R + 2A fc ) h T »» + -h\ a [u - |) 



3 

2(4 - a) m / m R\ 2 - a „ 6a . . 

+ 2 (2a -3f> +^) aa ( D+ R\ h _ h (w - *r+ Ifa - i D + B + a* , h 



9a V 2 / V 72 " 18 3 

Note the factor 1/4 instead of the traditional 1/2. The additional, irrelevant factor of 
1/2 has been extracted in order to have an operator where the coefficient of □ acting 
on h T is one. Note that the CS term affects only the propagation of }\F . Similarly 
we decompose the ghost field as 

C, = V, + d,S, (3.11) 
where V^V^ = 0, and similarly for C^. This leads to 



^ghost — / (Pz\r~g 



\ '" ( □ + — j k - ^—^sn ( □ + , 4 n , R \ s 

3 J " 2 V 3(3-/3) 



(3.12) 



Changing ghost variables to (V^, 5) again produces a Jacobian which will be canceled 
by another Jacobians coming from a redefinition of S (see (14.2j) ). 
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4 The Cutoff 



In this section we define the cutoff and we set up the calculation of the traces in (13.61) 
by means of harmonic expansions on S 3 , viewed as Euclideanized de Sitter space. 
From here on we will work in the "diagonal" gauge in which a — h couplings vanish. 
This fixes 



P 



2a + 1 



Furthermore we make the field redefinitions 



R 

□ + j = 4 ; 



f □ + ^ ) rr -.a: VQS' = S 



(4.1] 



(4.2) 



whose Jacobian determinants cancel those coming from (I3.9P and (13. lip [T6l[T7]. Then 
the action (I3.10p becomes 



S {2) + S GH = ~Z k I d s xy^ h TT ^A 2 , u ^h T J + + c a aA a a + c h hA h h 



(4.3) 



where we have defined the operators 



A P a 



(□ - f + 2a) 5{:^ + ^(/% ) v AV /n(n-f) 



A 



in 



A r 



□ + ^r^R + 2aA , 

,_, 2 — a „ 6aA 

□ + - i? + 



A, 

and coefficients 



□ + 



4- a" ' 4- a ' 
i2 6A 



A - a A - a 



c x = -- 



2(4 - a) 



A- a 



a 9a 18 

Similarly, the ghost action (I3.12p in the diagonal gauge becomes 



c(2) 
° ghost 



d 3 xy/—g 



V^Ay^Vy + c s S f □ + t-^ — R I ,S' 



(4.4) 



(4.5) 



(4.6) 



V A - a 

where c$ = {a — 4)/3. It would be tempting to choose a = 1, which corresponds to 
the de Donder gauge and simplifies the inverse propagator, but we shall not do so 
in order to probe the a-dependence of our results. We shall come back to this issue 
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later. We note that the gauge a = 4 which has been used in [2] is not singular, as 
can be seen from (13.1 Op . However, in this gauge the <J 2 ,h 2 , S 2 terms do not contain 
□ and the method we shall use below to compute the beta functions will turn out 
to be not suitable to deal with this case. Therefore, in the rest of this paper we will 
assume that a^4. 

Now we are ready to discuss the choice of the cutoff. For each spin component we 
choose the cutoff to be a function of the corresponding operator given in (14. 4p . Then, 
defining the gauge fixed inverse propagator 

/A, \ 



o 



ciAi 



C a A a 



\ 



(4.7) 



c h A h ) 



we choose the cutoffs to have the following forms 

/ R k (A 2 ) 



\ 



(4- 



cii4(Ai) 

c a R k (A a ) 

\ c h R k (A h ) J 

where R k (z) is a real profile function with the properties described in Section 1. Then 
the denominator of the first term in (13.61) is 

( P k (A 2 ) \ 
ciPkiAt) 

\ C h P k (A h ) J 

where we have defined the function P k {z) = z + R k (z). The bare couplings appearing 

in the second variation of the action do not depend on k, so in the numerator of (13.61) 

we use the following formula 

d 



o + n k = z 



(4.9) 



dt 



(ZaR k (Ai)) = z Ci d t Rk(Ai 



(4.10) 



Then the overall factors Z and q cancel between numerator and denominator and 

(13.61) reduces to 
dT k 



te- 



eth 



- [Tt 2 W(A 2 ) + TW(A0 + Tr H/(A CT ) + Tr W(A h )\ 

- [T^WiAv) + Tr ^(A 5 )] , 



(4.11) 
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where we have defined the function W(z) = 

The traces are given by the sum of the functions over all the eigenvalues. Such 
sums can be performed after Euclidean continuation to a three sphere. This is de- 
scribed in Appendix B. See Section 5 for further comments. 

We must now specify the cutoff profile R k . In order to be able to perform the 
summations in closed form we will use the "optimized" cutoff Rk{z) = (k 2 — \z\)9{k 2 — 
\z\) [IB] . Note the absolute value, which is normally not used because the operators 
one normally considers are positive. Insofar as the role of the cutoff is to limit the 
sum to a finite number of modes, this definition will work also if the operator is 
negative. We will discuss this point at some length in the next section. We have 
d t R k (z) = 2k 2 9{k 2 - \z\) and for z < k 2 , P k {z) = k 2 , so that W(z) = 26(k 2 - |z|). 
Dividing numerator and denominator by k 2 , they are given by 

± n n 

+ ^m^(l-A-) + ^m^(l-^) 

n n 

-2 £ 0(1 - )£) - 2 £ m S n 6(1 - X s n ) . (4. 12) 

n n 

Here X^ = X^ /k 2 are the distinct dimensionless eigenvalues of the Euclideanized 
operator A i; and rrin their multiplicities. Using the results of Appendix B, the 
eigenvalues are given by 



R 



/i(e) 



3/2 



= -(n 2 + 2n + 2) - 2A ± - ( ^ ) ;/(;/ + 1 )(// + 2 ) . // > 2 
R 

A« = — (n 2 + 2n - 3 + 2a) - 2aA , n > 2 , 
6 

R ( , 6(2 - a) \ 6aA 

\l = - n 2 + 2n - - , n > 2 , 

6V 4-a J 4-a ' ~' 

,h R ( 9 6 A 6A . , 

and the respective multiplicities are 

= = n 2 + 2n — 3 , 

m n = m n = 2(n 2 + 2n) , 
m ' J n = m ^ = m s n = n 2 + 2n + l. (4.14) 
12 



We see that the effect of the cutoff function is simply to terminate each sum over 
eigenvalues at some maximal value n max . In the next section we discuss the way in 
which one can evaluate these finite sums, emphasizing the subtleties that arise in the 
spin two sector. 



5 Evaluation of the functional Traces 

The most familiar way of evaluating the asymptotic expansion of functional traces 
is the method of the heat kernel. However, in the presence of the /i-dependent 
terms in the trace over the spin 2 fields, the heat kernel coefficients are not known. 
Since on the other hand for a maximally symmetric background the eigenvalues are 
known explicitly, we will directly evaluate the expansion of the traces using the Euler- 
Maclaurin formula. In the present context, where the function to be traced and all 
its derivatives vanish at +00, it says that 

00 r°° 1 r r 

E F ^ = dx F(x) + -F(no) - F'Cno) - F"'(n ) + R (5.1) 

n=no n ° 

where B n are the Bernoulli numbers and R is a remainder. For the evaluation of 
the beta functions of interest, the terms with three or more derivatives of F give 
no contribution, so it is enough to include the term containing B2 = 1/6 and the 
remainder can be neglected. 

The evaluation of the sums can be done using algebraic manipulation software. 
There are some subtleties that need to be emphasized. For each type of field, the 
effect of the step function is that the integral extends up to a finite value of n. Thus, 
in any spin sector in (I4.12p we have to evaluate, schematically 



ni ax 



dxm(x)6(l — A(x)) = / dxm(x) , (5-2) 

J n 



where n max is a (real and positive) root of the equation 



X(x) -1 = 0. (5.3) 



5 In general n max is not an integer. One can think that the sums can be evaluated by terminating 
them at the largest integer n < n max . However, it is easy to check that his procedure applied, for 
example, to the evaluation of the known heat kernel coefficients will not give the right answer. 



13 



When the eigenvalues are quadratic in n, the equation (15.3j) has two roots. Since the 
cutoff must be positive, it is always clear which one to choose. For a < 4, which 
is the case for the gauges we use, the trace h has a negative definite kinetic term. 
This is the well known problem of the indefiniteness of the gravitational action. In 
the evaluation of the sums (I4.12p this sign has no effect. One can choose n max as for 
the other components and this does not lead to any pathology in the flow equation 
(see [15] for a discussion). 

Things are more complicated in the presence of the CS term. The eigenvalues of 
the TT field contain a term which is cubic in n, and this term occurs with opposite 
signs in the two chirality sectors. If we tried to define directly the one loop effective 
action, this would lead to convergence problems similar to those discussed in [T9] . 
For us, convergence should never be a problem because we are only interested in the 
fc-derivative of the one loop effective action, which is finite due to the cutoff term. 
Still, one has to deal with the fact that the + and — eigenvalues have different signs 
(at least for sufficiently large n) and grow at different rates. 

Since for A T± equation ( 15 .3p is cubic, one has to choose one among its three roots 

r - s + t ■ r + e iT/3 s - e~ in/3 t ; r + e~ in/3 s - e in/3 t , (5.4) 

where r, s, t depend on R = R/k 2 , A = A/k 2 and ft = \xjk. Superficially the first 
root seems to be real and the others complex, so one could try to define the cutoff 
using the first root. However, the expressions s and t actually involve cubic roots and 
there are regions in parameter space where these roots develop imaginary parts. In 
fact, there is no root which is real for all parameter values. 

One can understand this better by considering the plot of the functions A T± (x) 
for x > 0, see Figures 1,2 and 3. Consider first the case when ft > 0. The function 
X T+ (x) is monotonically increasing, and for this function equation ( 15.31) has a single 
real root. Thus there is no ambiguity for the positive chirality modes. In the case 
of the negative chirality modes, however, for small x the function \ T ~(x) initially 
grows with x, until the cubic term prevails; hereafter it decreases. When R and A 
can be neglected, which is the situation we are interested in, the maximum is equal to 
Aft 2 /27 . Thus if ft > a/27/4 equation (15.31) has two positive roots. It is clear that the 
smaller of the two has to be chosen, namely the one where the function is growing. 
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Figure 1: The eigenvalues (solid curve) and (dashed curve) as functions of n, for 
R = A = 0.01. Left panel: large jl regime (here jl = 3). Right panel: small jl regime (here 
jl = 0.3). 



For this reason we will call this an "ascending root cutoff". On the other hand if 
< jl < a/27/4, equation (15.31) has no positive root and the ascending root does not 
exist. However, we observe that the equation \ T ~(x) = —1 has a real positive root 
for any jl. As mentioned in Section 3, in the evaluation of the beta functions it is not 
important whether the modes have positive or negative eigenvalue. Therefore we can 
use this descending root to define the cutoff. We call this a "descending root cutoff" . 
When jl < 0, the discussion can be repeated interchanging the roles of A T+ and A T ~. 

Since the descending solution of A T ~ = — 1 always exists, while the ascending 
solution of A T ~ = 1 only esists if |/2| > a/27/4, one wonders why not use always the 
former. The reason is that the information one can get from these different cutoffs is 
complementary. For example if we restrict ourselves to jl > 0, when jl becomes very 
large, the descending solution of A^ _ = — 1 is much larger than the ascending solution 
of A^ + = 1. In fact the sum over negative chirality modes is divergent in the limit 
ji — > oo. The beta functions computed in this way will have fictitious singularities 
in this limit and it will not be possible to flow smoothly to the case when the CS 
term is absent. As we will see later, this also implies that the beta functions have 
fictitious singularities near G = 0, so it will not be possible to study the Gaussian FP 
in this cutoff scheme. On the other hand, the ascending solution has an imaginary 
part when \jl\ < a/27/4 and therefore cannot be used to study the small jl region. 
The conclusion is that in order to properly understand the behavior of the theory 
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Figure 2: The real (left panel) and imaginary (right panel) parts of the roots of the equation 
A^j + = 1, for R = A = 0.01, as functions of fi. The first root (thick dotted line) always 
has positive real part and is complex for —^27/4 < fx < 0; the second root (dashed line) 
always has negative real part and is complex for < fl < a/27/4; the third root (solid line) 
is complex for —^27/4 < fl < \/27/4. The solutions of the equation = 1 are obtained 
by the reflection fi — > —ft. 




Figure 3: The real (left panel) and imaginary (right panel) parts of the roots of the equation 
A^~ = —1, for R = A = 0.01, as functions of fi. The first root (thick dotted) always has 
positive real part and is complex for fX < 0; the second root (dashed) always has negative 
real part and is complex for fi > 0; the third root (solid line) is always complex. The 
solutions of the equation = — 1 are obtained by the reflection fl — > —fl. 
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over the whole range of values of \x one has to use both schemes: the "descending" 
cutoff for small jx and the "ascending" cutoff for large jl. 



6 The Beta Functions 



Applying the Euler-Maclaurin formula to each of the sums in (I4.12p and extracting 
from each of the integrals the appropriate powers of R, the r.h.s. of ( 14.121) can be 
written 



d t Y k = 



C R- 3 / 2 + C 2 R^ 2 + C 3/2 + ~F(n ) - ^F" 



2! 



(6-1) 



where the sum is over h T , £, a, h, C T , C. 

The contributions of each spin component to the C3/2 term are independent of the 
cutoff scheme. They are listed in the following table: 





n 


C3/2 


F(n ) 


F'(n ) 


h T 


2 


12 


20 


24 


e 


2 


-24 


32 


24 




2 


-18 


18 


12 


h 





2 
3 


2 


4 




1 


8 
3 


12 


16 


C 


1 


16 

3 


8 


8 



In the second column, the lower end of integration in ( 15.11) is shown. Summing 
all the contributions we conclude that the coefficient of the i?-independent term is 
exactly zero. More precisely, the sum C 3 / 2 + |F(n ) — ^-F'(n ) is zero separately for 
the trace-free part of h^ v (the sum of the first three lines) , for the trace part h and for 
the ghosts (the sum of the last two lines). These are the components of the fields that 
can be defined by purely algebraic conditions. The cancellation does not occur for 
components defined by differential constraints, such as the transversality condition. 
To some extent, the overall cancellation of the R- independent terms is expected. In 
the simpler setting of pure gravity without CS term, or any matter field coupled to 
gravity, the sums can be evaluated using the heat kernel expansion. On a manifold 
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without boundary the trace of the heat kernel contains only even powers of R, and 
in a 3-dimensional manifold the volume prefactor is proportional to R~ 3 ^ 2 , so that 
the expansion of dtTk contains only odd powers of R, and there is no /^-independent 
term. So the CS term will not be induced, if one starts without it. In Appendix D we 
show that the sums done with the Euler-Maclaurin method explained above exactly 
reproduce the heat kernel results for 3D gravity without CS term at one loop. The 
fact that all the contributions listed in the table above exactly cancel when properly 
summed is therefore a nontrivial check of our calculation. 

To obtain the rest of the beta functional dtTk there remain to sum the contribu- 
tions of type C and C 2 in (j6.ip ; the final result has the following structure: 

V(S 3 ) 



d t T 



•a k 



16tt 



k*A(A,jl) + kB(A,ji)R + 0(R 2 ) , (6.2 



where we have inserted powers of k such that the A- and /^-coefficients are dimen- 
sionless. The volume of S 3 with radius £ is V(S 3 ) = 2tt 2 £ 3 with £ = \J~^- 

Equation ( 16. ip is an expansion in R, whose coefficients are functions of A, fi and 
k 2 . For reasons that will become apparent below, we shall restrict ourselves to the 
parameter region where A and R are of the same order. Therefore we shall also 
expand the coefficients in (16. ip in powers of A, namely in Cq we keep at most terms 
linear in A while in C 2 we only keep the A-independent terms. We will give explicit 
expressions for A and B later. 

Evaluating the (Euclidean version of the) renormalized TMG action ( 13. ip on S 3 
background, it can be written in the form 

r * - ^ 3 > {^g-iL r+ wb^ 3/2 + °< fl! ») • < 6 - 3 > 

where we have used that the integral of the CS term on S* 3 is given by J tr(udu + 
|a; 3 ) = 32n 2 . The couplings A, G, \i are now renormalized couplings evaluated at 
scale k. Rescaling the coupling constants as 

G = Gk~ l , A = Ak 2 , n = fik, (6.4) 

so as to make them dimensionless, and comparing the t-derivative of (I6.3P with ( I6.2p . 
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we obtain: 




0/A - -77- A I -i^ + A, (6.5) 



(6.6) 



8ttG I G y 8ttG 16tt ' 

aG 1 b 




16vrG 2 16ttG 16tt 

= 0, (6.7) 



12y/6irp,G 

The last equation results from the fact that the terms of order i? 3//2 in (16. 2 j) cancel, 
and it implies that the dimensionless combination z/ = G/i = G/i does not run. From 
the other two equations one obtains the one loop beta functions of G and A: 

d t G = G + 5(/a)G 2 , 

1 X 



d t A = -2A + -G (a(/2,A) +2S(/i)A) . (6.8) 



These equations have exactly the same form as in pure gravity with cosmological 
constant, except that the coefficients A and B are now /i-dependent. In order to 
solve for the RG flow, we use that v does not run and substitute jx = vjG. This 
yields two ordinary first order differential equations for Git) and A(t), depending on 
the fixed constant value of the external parameter v. Rather surprisingly, despite the 
very different functional form, the resulting flow is numerically quite similar to that 
of pure gravity with cosmological constant. Unlike in pure gravity, however, here we 
cannot give the solution in closed form. We will now describe the results for different 
cutoff schemes. 



7 Results 

7.1 Ascending root cutoff 

For \p\> we define the cut on the spin two modes as the smallest root of the 

equation A T± = 1. In the case of positive p,, the cutoff on A T+ corresponds to the 
thick dotted line in Fig 2, while the cutoff on A T_ is obtained by reflecting the solid 
line. (Note that the two roots are different, so that the number of modes of positive 
and negative chirality is different.) Calculating the sums yields beta functions that 
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are real and well defined for < G < a/4/27za For negative jl the roles of A T+ and 
A T ~ are interchanged. Calculating the sums yields beta functions that are real and 
well defined for — a/4/27z/ < G < 0. The two calculations match smoothly along the 
line G = 0, so one can put them together to obtain a RG flow on the whole region 
G 2 < Au/27. 

With this prescription we calculate the coefficients A and B. These arise as 
complicated functions involving cubic roots, but after some manipulations can be 
reduced to the following relatively simple form: 

~ 16 9(2 v / 3cos2^- yft cos AO + 8(cos#) 3 sin 9) 
A(A, n) — — - — I —To (7.1) 

8(3 + lla-2a 2 )^ + 48(cos#- \/3sin#)^ 
7r(4 — a) -KsinQO 

4(l + a)(ll-2a) 2( v / 3 sin 9 - cos 6) + 22( v / 3 sin 59 + cos 56)) 



B{p) = 



7r(4 — a) 37rsinQ9 
where we have introduced the angle 



1 Au 2 
= -arctaay-^--l. (7.2) 

The beta functions admit a Taylor expansion around A = G = 



~ - 2(94-4a 2 )G 2 95G 4 2233G 6 / - 7 x 
^ G " G " 37r(4-a) "fa?"&7 + l G J ' (? 

4(14- 27a + 4a 2 )GA G 3 (42 + 115A) 11G 5 (78 + 343A) ^ 6 



3tt(4 - a) ' 6ttz/ 2 32ttz/ 4 

= 4(14-27 n + to»)6A + G^(42 , ;-' +0 (g 

3vr(4-a) 6ttz/ 2 32ttz/ 4 V 

We observe, as a nontrivial check, that in the limit Ji — >■ oo the beta functions 
agree with the result for gravity without CS term, which are calculated in Appendix 
D using a different method. 

The flow is shown, in the case v — 5, in Figure 4. We note the main features, 
which hold for any value of v > 0. There is a FP, called the Gaussian FP, in the origin, 
which is seen to be attractive in the A direction and repulsive in the G direction. To 
make this statement more quantitative, one has to study the linearized equation 

k d ® i ~ § * ) =M iS (g j -~g jit ). (7.4) 
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-1.0 -0.5 0.0 0.5 1.0 -0.10 -0.05 0.00 0.05 0.10 



Figure 4: The flow in the A-G plane for a = 0, v = 5, in the ascending root cutoff scheme. 
Right panel: enlargement of the region around the origin, showing the Gaussian FP. The 
beta functions become singular at |G| = 1.9245, outside the domain of the picture, but this 
singularity is an artifact of the scheme. 



where 

M «=|L> ft = fc f- < 7 - 5 > 

Here gi are the couplings and <ji* their values at the FP. At the Gaussian FP = 
one finds 

The eigenvalues of this matrix ( "scaling exponents" ) correspond to the canonical mass 
dimensions of A and G, as expected. An important fact (which had been previously 
observed in the absence of CS term in [9]) is that the eigenvectors coincide with 
the coordinate axes. This is not the case in higher dimensional gravity, where the 
eigenvector with eigenvalue equal to the canonical dimension of G has a novanishing 
component along A. 

In addition to the Gaussian FP there is a nontrivial FP at A* = 0.000490471 and 
G* = 0.200016. It is seen to be UV attractive in both directions. The eigenvalues of 
the corresponding stability matrix are found to be equal to —2.29401 and —1.00515. 
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Figure 5: Position of the FP (left panel) and eigenvalues of the stability matrix (right 
panel) for the nontrivial FP with a = 0, 1 < v < 40, in the ascending root cutoff scheme. 
Note that for this range of v the singularity is always above the FP. In the left panel, v 
grows from right to left. Note that A* > in this scheme. See Section 7.1 for a discussion. 
For large is, G* tends to 0.2005 and the eigenvalues tend to —1 and —2.298. 



The position of the FP and the eigenvalues of the stability matrix are given, for 
other values of z/, in Figure 5. Note that A is always positive but very small. This is 
due to the absence of a term of order G 2 in the expansion of dtA in (17. 3ft . Actually 
if one plots the contours of the beta functions for A one may see this as an effect 



of the deformation of the flow due to the presence of the boundary at G = a/4/27 
(which occurs at G = 1.924 in Figure 4). This can probably be regarded as a scheme 
artifact. 

Finally we observe that in the limit jl — > oo the beta functions reduce to 

Sty 4 — a 

A 14 _ 27<y 4- An 2 ~ ~ 
d t A = -2A-— Z ' a + ^ a AG. (7.7) 
3n A — a 

which agree with the result for pure gravity without CS term. 
7.2 Descending root cutoff 

In this scheme the cut on the positive and negative chirality spin two modes is defined 
by two different equations, as follows: 
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fx > 


fx < 


A T + = 1 
A T " = -1 


A T + = -1 



Using the criteria described in the previous section, in the case of positive fx, the 
cutoff on A T+ corresponds once again to the thick dotted line in Fig 2, while the 
cutoff on A T ~ is the thick dotted line in Fig 3. For negative fx the roles of A T+ 
and A T ~ are interchanged. Again the two roots are different, so that the number of 
modes of positive and negative chirality is different. In fact in this case the behavior 
is drastically different, since the cut on the descending mode grows linearly with fx. 
This will give a singularity in the beta functions for fx — > oo, and therefore, when we 
make the substitution fx — v/G, for G — >■ 0. Thus this scheme is really useful only 
for sufficiently small fx. 



With this prescription we find: 

2 cosh 2?7 



A(/2,A) 



16 y/3 

3n 7r 



cosh 3i] 

(3 + 11a -2a 2 ) ~ , 8^ 



2 cosh 20 + 1 



cosh 3ip 



A + 



B(fx) 



+ 



7r(4 - a) 

4(11 + 9a -2a 2 ) _ 
3tt(4 - a) 

2 / 8cosh2r7 



7T 



I 



+ 



I 



20^3 



2 cosh 77 + cosh 3i] 2 sinh ip — sinh 3ip 

cosh 2i] + cosh 2ip 
37r cosh 3i] 

8 cosh 2-0 + 1 



(7.8) 



A 



+ 



3 \Z3tt V cos h 3^ + 2 cosh 77 sinh 3ip — 2 sinh ip 
where we have defined 



77 = - arctanh \ / ! 



4/i 2 
~27 



, 1 1 /- ^ 2 

ip = - arccoth d 1 + — . 



(7.9) 



A representative flow is shown in Figure 6 for v = 0.1. Superficially this may look 
similar to Figure 4, but there are some important differences. First and foremost, 
there is no Gaussian FP, as is clear from the enlargement of the area around A = 
G — 0. The FP is wiped out by the singularity at G — 0, where the beta function of 
A blows up. Note that the flow lines run in opposite directions on the two sides of 
this line. All the flow lines which, for positive G, seem to come in from plus infinity 
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-1.0 -0.5 0.0 0.5 1.0 -0.10 -0.05 0.00 0.05 0.10 



Figure 6: The flow in the A-G plane for a = 0, v = 0.1, in the descending root cutoff 
scheme. Right panel: enlargement of the region around the origin, showing that there is no 
Gaussian FP. The beta functions diverge on the A axis. 

actually originate from minus infinity, and follow closely the G axis before bending 
over. Similarly for G < all the flow lines that seem to come in from minus infinity 
actually come from plus infinity and run close to the G = line. There is thus an 
infinite accumulation of flow lines along the A axis. Another difference with Figure 
4 is that the flow lines near the G axis are tilted in the opposite direction, and by a 
much larger amount. 

For v = 0.1 the nontrivial FP occurs at A* = -0.0565337 = 0.30036 and the 
eigenvalues of the stability matrix are —2.76746, —0.781453. The position of the FP 
and the eigenvalues, for 10 -6 < v < 0.5 are shown in Figure 7. The FP occurs at 
positive A for v > 0.18 and negative A for v < 0.18. Given that this cutoff scheme is 
more dependable for small jl (and hence, at fixed G, for small z/), it is again possible 
that the positive values of A are a scheme artifact. 

Just as the beta functions in the ascending root scheme tend to a very simple form 
in the limit jl — > oo, the beta functions in the descending root scheme tend to a very 
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Figure 7: Position of the FP (left panel) and eigenvalues of the stability matrix (right 
panel) for the nontrivial FP with a = 0, 10~ 6 < v < 0.5, in the descending root cutoff 
scheme. In the left panel, v decreases from right to left. The cosmological constant changes 



sign for v = 0.18. The rightmost point [y = 0.5) has p ~ 3 > y27/4 and therefore is in 
the region where the scheme becomes unreliable. 



simple form in the limit p — > 0: 

= 6 _ ± l 1 + 9 a -2a^ 
37r 4 — a 

8,A = -2A-i-fl + i^±^ A)a. (7,0) 
37r \ 4 — a / 

These beta functions have a FP at A* = — 1/3, = ^-^^-2^) ' wm ch for a = is 
~ 0.8568. The eigenvalues of the stability matrix are —2.182 and —1 corresponding 
to the directions (1,0) and (-0.5499,0.8352). 



7.3 Spectrally balanced cutoff 

In the preceding sections we have described two calculations of the beta functions for 
TMG. Both schemes are "spectrally unbalanced" in the sense that the summations 
over the spin two fields in (14.121) contained a different number of positive and negative 
chirality modes. Consideration of the behavior of the roots in Figures 1, 2 and 3 show 
that the ascending scheme becomes balanced for large fl while the descending scheme 
becomes balanced for small p. These are the regimes where these schemes are most 
reliable. One thus wonders whether one can work in a scheme where the sums always 
contain equal numbers of positive and negative chirality modes, by construction. 
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Such a scheme can be constructed by tweaking the cutoff profile. First, we must 
allow different profiles for different spin components. The spin one and zero com- 
ponents will be treated as before. For the spin two components we allow the more 
general form Rk{z) = (q(z)k 2 — z)8(q(z)k 2 —z) where q(z) is a dimensionless function. 
For any choice of the function q(z) we have d t Rk{z) = 2k 2 6(q(z)k 2 — z) and for z < k 2 , 
Pk(z) = k 2 , so that W(z) = 29(q(z) — z/k 2 ). With this cutoff choice the traces are 
still finite sums of the multiplicity, up to a maximal value n max which is determined 
as a root of the equation 

\( x ) = q(X(x)) . (7.11) 

Thus the sums are still calculable by the same technique used previously. A choice 
of n m ax implicitly defines a choice of the function q. The sums on the spin two 
components in (14.121) are now replaced by 



n 



max 



EE^WCA^) - A^) = 2 £ ml ■ (7.12) 



± n=2 n=2 



If Jx > 0, we choose n max to be the unique positive root of the equation = 1; 
then we cut both sums at this value. In this case the function q is q + = 1 for 
the positive chirality modes, while q~ is determined by the condition g~(A^~(x)) = 
A^ + (x) I (x) . Likewise, if jl < we cut at the root of = 1. 
With this prescription we get the following results: 

^• A > - -l + ^kf (7 - 13) 

8(3 + lla-2a 2 U 16^ ~ 
-A + — ^ , 



7r(4 — a) 7r(cosh 3r] + 2 cosh 77) 

4(11 + 9a -2a 2 ) 8^/ 8 + llcosh2r? 
37r(4 — a) 9n \ cosh 3rj + 2 cosh rj 



where r\ is given by (17. 9p . Thus, we observe that the apparent poles at a = 0, which 
correspond to 77 = 0, are actually absent as can be readily deduced from (I7.13p . and 
the coefficients A, B are real. Note also that 77 = i9 where 9 is the angle that we 
have encountered in Section 7.1 (see (17. 2\\ . Furthermore cosh 37/ = 3\/3/|Al> and the 
results depend on the absolute value of jl. 
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In the limit — > oo we get 

8(11 + to - 2^ 4(47-2a-) 

7r(4-a) 3(4 -a) v ; 

This limit corresponds to neglecting the CS term, and it is reassuring that we find 
agreement with the result of the heat kernel calculation in Appendix D. 

The Taylor expansion of the beta functions in powers of G around the Gaussian 
FP is 

9t G = G- 4(4 3 7 - 2 " 2 f + fg-^ + (ffl), (7,5, 
37r(4 — a) Sttv bTiu 2 V / 

d t K = 2A 4(14 ~ 27a + W)dK 4G 2 (3 + 5A) | G 3 (42 + 115A) | Q ( & 

37r(4 — a) 3ttu Qtxu 2 V 

The ^-independent terms coincide with the result for pure gravity, as described in 
Appendix D. 

In the limit p, — > the beta function of A has the same expression as in (I7.10p . 
but the beta function of G becomes singular. Still, the position of the FP seems 
to approximate closely the one that was found in the descending root cutoff. For 
v = 10 -6 we find A* = —0.315, = 0.830. This should not come as a surprise, since 
in this limit the two positive roots of the equations A T+ = 1 and A T ~ = — 1 become 
equal, so the descending root cutoff becomes spectrally balanced. 



8 Discussion 



8.1 Summary 

We begin by summarizing the main results. We have calculated the beta functions 
of A, G and fi and we have studied the corresponding RG flows. It is generally the 
case that the beta functions of dimensionful couplings (by which we mean also the 
beta functions of the corresponding dimensionless ratios that we denoted by a tilde) 
are scheme dependent even to leading order. This is in contrast to the more familiar 
case of the dimensionless couplings, whose beta functions are scheme independent 
to leading order. It is therefore not too surprising that the three schemes we have 
considered give different pictures. Out of these different pictures one can however 
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Figure 8: Left panel: The position of the FP in the A-G plane, varying v from 0.002 (upper 
left) to 1000 (lower right). The point with coordinates (0,0.2005) is the limit v —> oo. Right 
panel: The eigenvalues of M as functions of v. For v = 0.002 they are —1 and —2.298 while 
for v = 1000 they are -0.969 and -2.238. 

extract some common features that are presumably scheme- independent. We sum- 
marize them here. 

The first main result is that the dimensionless combination v = fiG, which is the 
coefficient of the CS term, does not run. We will comment further on this in the 
subsection on topological properties. Due to this fact, the remaining beta functions 
describe a flow in the A-G plane, whose properties depend on the value of the fixed 
constant v. This two dimensional flow is governed by two FPs: a Gaussian FP in 
the origin and a non-Gaussian FP at positive G. The Gaussian FP is UV attractive 
in the A direction and UV repulsive in the G direction, with scaling exponents equal 
to the canonical dimensions of A and G. This FP is not seen in the descending root 
cutoff scheme, but this is clearly an artifact of the unphysical singularity of the flow 
for large \x (hence small G). For the same reason one should not trust the properties 
of the Gaussian FP that reappears in that scheme in the limit v — > 0. 

There is then a nontrivial FP which always occurs for positive G. Whether this 
FP has positive or negative cosmological constant seems to be somewhat unclear. The 
ascending root cutoff, which is reliable for large u, and the spectrally balanced cutoff 
both indicate that for large u the FP tends towards the FP that is known to occur 
in the absence of CS term, which has A* = 0. The descending root cutoff, which is 
reliable for small v, and the spectrally balanced cutoff both indicate that for small v 
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the FP tends towards the values A* w —0.33, G* ~ 0.85. We conclude that for small 
v the cosmological constant is certainly negative. The fact that for v > 0.18 in the 
descending cutoff scheme A* becomes positive is probably an unphysical effect of the 
singularity occurring at G = 0, as Figure 7 suggests. Whether for intermediate values 
of v A* could be positive is an issue that we leave for future more refined analyses. 
On balance, we think that Figure 8 probably gives the most reliable picture. 

The final universal feature of the flow is that the nontrivial FP is UV attractive 
in both directions, and the existence of a trajectory that connects the nontrivial FP 
in the UV to the Gaussian FP in the IR. The scaling exponents of the nontrivial FP 
depend on v but both for very large and very small v they tend to similar values 
which are close to -1 and -2.3 (exactly -1, for large v). 

8.2 Small G expansion 

We have confined ourselves to the one loop approximation, which is the lowest order 
in an expansion in G. From this point of view one could be surprised by the fact that 
the beta functions seem to contain arbitrarily high powers of G (see equations (17. 3p . 
( 17.151) ). But note that the original beta functions of A, G and fx are indeed at most 
quadratic in G: It is only when we use the constancy of v to solve for jl as a function 
of G that the beta functions become nonpolynomial in G. Since the beta functions 
will anyway receive higher order corrections, one may want to truncate the expansion 
in ( I7.15P to terms at most quadratic in G. 

In the case of the ascending root cutoff there are no terms of order G 2 beyond 
those that are already present in the absence of the CS term, so to leading order in 
that scheme the FP is the same as in the absence of CS term. In the case of the 
descending root cutoff one cannot even pose the question, since the beta functions 
are singular in the origin. In the case of the spectrally balanced cutoff one obtains 

^ = «_1(£_^, (8 ,) 

37r(4 — a) 

~ ~ 4(14- 27a + 4a 2 )GA 4G 2 (3 + 5A) 

OfI\ — — 2J\ — - — . 

37r(4 — a) 3nu 

In this case the beta function for G is the same as for pure gravity without CS term, 
but the beta function of A receives a ^-dependent correction. In this approximation 
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the location of the FP is at 



G„ = — ; A* = — . (8.2) 

47 ' 5tt + 423z/ V ; 

and the eigenvalues of the linearized flow equation are —1 and — ^ — ■ 

In any case we believe that the vanishing of the beta function of v will be true also 
when higher order corrections are taken into account. This provides some justification 
for keeping also the higher order terms in the beta functions. 

8.3 The chiral point 

One of the motivations for this work was to determine whether the chiral point defined 
by the condition /x = ±a/|A| has any special properties under RG. We can recast this 
issue in a slightly more general way as follows. With the dimensionful couplings A, 
G and \x one can construct the dimensionless combinations (|3.2p . The more general 
question is whether any of these combinations is invariant under the flow. We have 
seen that this is the case for v. The conditions = const, r = const determine 
two dimensional subspaces in the space of couplings. The intersection of these sub- 
spaces with the A-G plane defined congruences of curves. The lines of constant r are 
hyperbolas with equation 

G 2 = I (8.3) 
A V ; 

Since v is RG invariant, using to the relation 

we see that also is constant on these hyperbolas. It is easy to see that the flow 
lines will generally not preserve 0: they will intersect the lines of constant 0. If we 
choose the initial conditions for the flow to lie on one hyperbola, then the flow will 
move away from it. In this sense, the condition = 1 is not a RG invariant. 

One can still ask whether there exist some FP that satisfies the chirality condition. 
For each value of u, the corresponding FP in the A-G plane determines also a FP value 
for jl; with these FP values, one can compute the dimensionless combinations is, r and 
0. With the spectrally balanced cutoff, there exists a unique value is crit ps 0.0956016, 
for which the FP lies exactly on the corresponding hyperbola. For this choice the 
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chirality condition is UV attractive, but even in this case the flow does not preserve 
the hyperbola. 

The only point in the A-G plane where <fi becomes approximately RG invariant is 
near the Gaussian FP, where the dimensionful A, G and /x become constant. In terms 
of the dimensionless A and G this means that in their beta functions one only keeps 
the first term on the r.h.s. of (17.71) . Then the flow lines approximate the hyperbolas 
of constant r (see right panel in Figure 4). The approximation becomes better when 
t is small. If we want to have = 1 approximately constant this can be achieved for 
v = ^ ^ 0. 

8.4 On-shell properties and the sign of A 

Given that the effective action should be gauge independent on-shell, i.e. for R = 6A, 
also its t-derivative should have this property. From equation (16. 2p . this implies that 
A + 6-BA has to be a- independent. Indeed one finds that 

A + 6-BA = -1^— (1 + 3A) + /i-dependent terms (8.5) 

so the a-dependence cancels. Another way to see this is to observe that on-shell the 
effective action is proportional to 1/y/r = 1/ (Gy/X), up to terms of order yX. Thus 
the beta function for this particular dimensionless combination of couplings should 
also be gauge independent. Indeed, using (16. 8 j) one finds that this beta function is 
proportional to A + 6-BA, and therefore a- independent. In view of (18. 4p . the same 
is true for the beta function of = fi/y/\A\. Note however that the beta functions 
of t and <p are not functions of r and/or alone but depend on A and G separately. 
Therefore, to calculate their flow one needs the beta functions of A and G, which 
have to be computed with off-shell backgrounds. 

It is perhaps appropriate to add some further comments on the significance of 
the on-shell condition. One may worry that off-shell backgrounds will give rise to 
corrections proportional to the equations of motion. We show in Appendix C that 
this is not the case. One may also worry that since the sums (14. 12j) are performed 
on a three sphere, which has positive curvature, the resulting beta functions are only 
correct in the domain A > 0. In particular, could this not affect the existence of the 
FP, which actually occurs (in the descending and balanced schemes) for negative A? 
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Would one not have to perform the calculation on an AdS background? It would in- 
deed be desirable to perform such a calculation, but since the Euclidean continuation 
of AdS spacetime is a noncompact hyperboloid, the harmonic analysis would be much 
trickier. Fortunately there are good reasons to believe that in the UV limit such a 
calculation would yield exactly the same beta functions that we have derived here. 
The reason is that any positive definite metric looks the same at sufficiently short 
distance scales, so the beta functions would be the same as long as k 2 3> R. Further- 
more we observe that for large v the actual values of the cosmological constant, in 
units of k, can be made arbitrarily small and negative, while the background can be 
chosen to have R arbitrarily small and positive, so there exist FPs which are as close 
to the mass shell as one wishes. 

8.5 Topological issues 

The vanishing of the beta function of the coefficient of the CS term calls for a topologi- 
cal interpretation. The analogous phenomenon in gauge theories has been investigated 
in [291130]. It was found that in spite of the finiteness of the theory, the coefficient 
of the CS term receives a finite renormalization proportional to the Casimir of the 
gauge group in the adjoint representation. This however is an infrared effect [31]; in 
the context of the flow of the effective average action, which we used here, one can see 
that the coefficient of the CS term is constant for all finite k except for a finite jump 
in the limit A; — ?- [32J. Since our beta functions are only valid in the case k R, 
we do not see such effect and we cannot exclude that such a finite renormalization 
happens. 

A related point is the following. It has been argued in [33] that as long as one 
works on a sphere, it is not necessary to impose quantization of u, and since we use 
precisely a sphere as a background, it is not clear that the vanishing of the beta 
function of v has a topological explanation. The answer is that the beta functions we 
obtain are actually insensitive to the topology of the manifold and therefore should 
respect the quantization of v in those cases when the topology demands it. 
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8.6 Asymptotic safety and higher derivative terms 

In four dimensional gravity, a nontrivial FP has been found in the leading order of 
the 2 + e expansion [20] and of a 1/N expansion [21] , in perturbation theory [22|[23] 
and in a variety of truncations of an exact RG equation [21H26] also in the presence 
of matter [27]. The nontrivial FP we have found for each value of v in TMG is 
similar to the FP that is found in the Einstein-Hilbert truncation of gravity also 
in higher dimensions, except for the fact that the cosmological constant turns out 
to be negative. There are however two features that make our result particularly 
interesting. 

The first is that for large fx the FP value of G is securely within the perturbative 
domain. There is therefore every reason to believe that the one loop calculation we 
have performed correctly captures the main features of the Wilsonian flow of this 
theory. The analogous result in pure three dimensional Einstein Hilbert gravity is 
physically less appealing because in the absence of CS term the theory does not have 
propagating degrees of freedom. 

The second has to do with the consistency of the truncation. In principle the 
effective average action will contain also higher derivative terms, which in our com- 
putation have been neglected. However, in three dimensions the Riemann tensor is 
expressed entirely as a function of the Ricci tensor, and consequently it is possible 
to use field redefinitions to eliminate higher derivative terms by means of field redef- 
initions, order by order in perturbation theory [28]. In three dimensions the higher 
derivative operators are technically redundant and therefore the truncation consid- 
ered here is consistent. At least in perturbation theory, it is not necessary to consider 
an infinite set of operators. However, the redefinitions used to eliminate the higher 
derivative terms change the cosmological constant and Newton's constant. Therefore, 
the beta functions of A and G will receive corrections which have not been considered 
here. 
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A Useful lemmas 



In implementing various cutoff schemes, it is useful to work out separately the de- 
composition of various terms occurring in the quadratic action + Sqf given in 
( 13. 7p . The decomposition of the terms that come form the Einstein-Hilbert action 
with a cosmological term follows: 
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The decomposition of the terms arising in the gravitational CS terms gives 
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B Euclidean continuation and harmonic 
expansions on S* 3 



The issue of defining the Euclidean continuation for a theory of gravity is quite subtle 
and we will not try to address it here in generality. In general one would have to 
think of the real (three dimensional) Minkowskian spacetime as being a section of a 
three dimensional complex manifold, and find another real section where the metric is 
positive definite, and where the action is real and possibly bounded from below. The 
Wick rotation procedure which is employed in perturbative path integral formulation 
of non-gravitational field theories runs into severe problems in gravity in which the 
notion of time is affected by difleomorphisms. Nonetheless, here we will adapt the 
analog of the standard field theoretic Wick rotation, since we are interested in one-loop 
beta functions, which are calculable from the gauge fixed quadratic part of the action 
when expanded about maximally symmetric backgrounds characterized in (13. 8ft . The 
nature of the Wick rotation still depends on the choice of background metric. For 
our purposes, it is convenient to work with the dS$ metric, which can be represented 
as 

d s 2 = -dp 2 + cosh 2 p(d6 2 + sin 2 6d(f) 2 ) . (B. 1) 

Upon Wick rotation p — > —ip, this turns into a metric on S 3 given by 

ds 2 = dp 2 + cos 2 p{dQ 2 + sin 2 6d<p 2 ) . (B.2) 

Thus in the path integral e lS ^ becomes e~ SE ^ 9E \ where 

S(g) = i J dpd9d(j)y^7j C\g) -»■ -J dpdBd^^~ E C[g E ] = -S E (g E ) . (B.3) 

With this procedure the Euclidean action differs from (13. ip just by an overall sign. 

To perform harmonic analysis on S 3 , we decompose all the fields that carry ir- 
reducible representations of the stability group SO (3) C SO(4), where £0(4) is the 
isometry group of S 3 . These fields admit harmonic expansion in terms of all the 
representation functions of 5*0(4) that contain the given SO (3) representation (see, 
for example, [34]). Thus, the divergence-free and trace-free symmetric tensor h? u , the 
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divergence-free vector £ and a generic scalar <p are expanded on 5 3 as □ 

oo 

W = E(^ )y & 2) w+^" a)y &" a) (^)) . 

n=2 

oo 

n=l 

oo 

m = T,^'°K q ,0) ( x ) . ( B - 4 ) 

n=0 

where g labels a unitary irreducible representation (UIR) of 50(4), the functions 
5ju/, g , 5^t,g and Y 0q are respectively the /iz/'th, /zth and Oth row and qth column 
of the Wigner functions for 50(4) UIRs labeled by the highest weight (n,s), with 
n > \s\. In expanding a field in spin j representation of 50(3) only the 50(4) 
UIRs that contain the given 50(3) representation occur, and hence the restriction 
n > j — \s\- Moreover, the 50(3) representation is contained only once. The 
expansion coefficients h^ n ' ±2 ^ are constants, and x labels the coordinates on 5 3 . The 
Wigner functions occurring in these expansions have the same symmetry, trace and 
divergence properties as the fields that the expanded fields. Note the presence of 
the zero modes in the form of the single constant mode Y^°°\ the six Killing vectors 
Yp 1 ^ 1 ' in the adjoint representation of 50(4) and the four conformal Killing vectors 
d^Y {10 \ These obey 

9fl Y (°°) = , V { , V, } F ( 10 ) = , V^ 1 ' ±X) = , (B.5) 

where the notation {/z, v\ means the symmetric and traceless part. The Killing vectors 
and the conformal Killing vectors together generate the conformal group S0(4, 1) that 
acts on 5 3 . 

Next, we make use of the fact that the group theoretical considerations yield the 

6 We have absorbed the normalization factor \J ^j+i) occurring in the expansion of a spin j held 
into the definition of the expansion coefficients [53], where d(n,s) is the dimension of an 50(4) 
unitary irreducible representation with highest weight (n, s). 
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following equations 



~OY^f\x) = f[(n+l) a -3]y^f)(x), (B.6) 
-UY^\x) = j[(n+l)*-2]Y<***\x), (B.7) 

v^tf^W = ±^v // f (ri + 1)£V ^' ±1)(:r) ' (R8) 

-^(x) = f[(^+l) 2 -l]Vo ( ; 0) (x), (B.9) 

V^'.r^) = ±«y|(n + l)f^r^(x). (B.10) 

with multiplicities as given in (14.14jl . In deriving these formulae, we use the well 
known formulae for the Casimir eigenvalues of the second-order Casimir operator 
Ci{n, s) for an SO (4) representation labeled by (n, s), and it dimension d(n, s) given 
by 

C 2 {n,s) = n(n + 2) + s 2 , d{n, s) = (n - l) 2 - s 2 . (B.ll) 
Furthermore, we use the formula 

- OY^ s \x) = | [C 2 (n, s) - C 2 (i)] Y^ s \x) , (B.12) 

where Y^.'^{x) is the Wigner function for the SO (4) UIR in the (n, s) representation, 
viewed as a matrix, restricted in its row to the lowest spin j representation it contains, 
and C2O) = j(j + 1) is the eigenvalue of the second-order Casimir operator for 50(3) 
in spin j representation. 



C The off shell one loop beta functional 

In the evaluation of the one loop effective action using the background field method 
one usually assumes that the background is a solution of the classical equations of 
motion. On the other hand in the derivation of the beta functions for gravity we have 
to evaluate d t Tk for off-shell backgrounds. This is necessary in order to disentangle 
the beta functions of G and A. One may worry that for off-shell backgrounds the 
simple equation (12 .7p is no longer correct. We show here that this is not the case. 
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In order to do this, we review the derivation of that equation within the back- 
ground field method. We consider a classical action S(4>) and expand = (p + r/ where 
the background is arbitrary. We Taylor expand 



5S 



1 



where A^ = |g 
now be written as 



5(0) = 5(0) + J — _ r] + -jj rjA(<P)ri + ... (C.l) 
_. The k- dependent generating functional analogous to (12. 3 p can 



e" Wfc N = J D V exp j-5[0 + 77] - AS k [rj\ - J jr/| , 



(C.2) 



such that 

5W k 



V=(v) = -jf- (C3) 

The background-dependent effective average action is the fc-dependent modified Leg- 
endre transform 

T$, ¥>} = W k [j] -Jjip- AS k [<p] , (C.4) 



where j has to be expressed as a function of (p, solving f )0.3[) . From here one finds 
that 

ST k . 5 2 T k (5*W k \- 1 

-jj— = -j - Rk<p ; inr = ~ ~ Rk ■ c - 5 

dip 0(fd(p \ ojoj J 

Let us evaluate these functionals at one loop. We have 

wf } [0, j] = 5(0) + ^Tr log (A* + R k (At)) -\ j J j{Af + RtiA^j . (C.6) 

From here solving equation (IC.3I) we get 

j = -(A$ + R k (A$))<p. (C.7) 

Therefore the one loop average effective action is given by 

T k 1 \^ v ) = S^) + ^T\og(A^ + R k (A^) + 1 - J j tpA^. (C.8) 

The main point to observe here is that the terms of the form j (pR k <p have exactly 
cancelled. Therefore, the only dependence on k is in the trace term, and the beta 
functional (12.71) is unaffected. In a gauge theory, the presence of the gauge fixing 
terms does not change this conclusion. 
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D Gravity without CS term 



Here we derive the beta functions of A and G using heat kernel techniques. This has 
been described elsewhere in a variety of gauges and cutoff schemes, but not using 
the cutoff scheme described in Section 3. For the sake of comparison with the beta 
functions given in Section 5 we give here this calculation. 

In the absence of CS term, and in diagonal gauge, all the operators appearing in the 
functional traces ( 14. lip are minimal Laplace-type operators of the form A = — V 2 1+E 
where E is a linear map acting on the spacetime and internal indices of the fields. In 
our applications it will have the form E = (q\R + q 2 A) 1+ where 1 is the identity in 
the space of the fields and q±, q 2 are real numbers. 

The trace of a function W of the operator A can be written as 

POD 

TrW(A) = W(\i) = / dsTiK(s)W(s) (D.l) 

Jo 

where Aj are the eigenvalues of A, W is the Laplace anti-transform of W(s) and 
TrK(s) = e~ sAi is the trace of the heat kernel of A. The UV behavior of the 
theory is governed by the lower end of the integration over s. We use the asymptotic 
expansion for s — > 0: 

Tr ( e " A ) = [A>(A)s-* + B 2 (A)s-l +1 + ... + B d (A) + B d+2 (A)s + ... 

(D.2) 

where B n = J d d Xy/gtrh n and b n are linear combinations of curvature tensors and 



their covariant derivatives containing n derivatives of the metric. Then (ID.lj) becomes 



W(A) = [Q f (W)B (A) + Q hl {W)B 2 {A) + ... 

+ Q (W)B d (A) + Q^(W)B d+2 (A) + ...], (D.3) 

where 

poo i poo 

Q n (W) = / dss- n W(s) = — - / dzz n - l W(z) . (D.4) 
Jo T ( n ) Jo 

The last equality only holds for n > 0. 

For the function W = of interest in the evaluation of (14.111) the integrals are 

Q 3/2 (W) = -^k 3 ; Qy 2 {W) = —7=k (D.5) 
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The only remaining ingredient that is needed for the evaluation of (16.21) are the 
coefficients Bo and B2 of the operators A given in (14.41) . These can be evaluated using 
standard methods (see for example Appendix B in and are listed in the table. 

The result is 
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7r(4-a) ' 37r(4-a) 

These results agree, for a — 1, with those reported in Section IV D of [llj. As 
discussed there, the values of A and B do depend on the cutoff scheme. This is 
normal for the beta functions of dimensionful couplings. It is reassuring that this 
scheme dependence does not change the picture qualitatively. Furthermore, there are 
certain aspects of the flow that are scheme independent. The first one is the beta 
function of the dimensionless coupling v, which is zero independently on the scheme. 
Another one is the beta function of the dimensionless combination r: it is equal to 

36AG 3 
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TV 



(D.8) 



in both cases. This accords with the gauge independence of the same quantity, as 
discussed in Section 6. Higher terms of the expansion in A will probably not be 
likewise universal. 

We conclude by giving the properties of the nontrivial FP in this scheme. It occurs 
at A = 0, G = 0.2005; its scaling exponents are —1 and —2.298, with eigenvectors 
aligned with the coordinate axes. 
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